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Abstract. Each conservation law of a given partial differential equation is determined (up to equiva- 
lence) by a function known as the characteristic. This function is used to find conservation laws, to prove 
equivalence between conservation laws, and to prove the converse of Noether's Theorem. Transferring 
these results to difference equations is nontrivial, largely because difference operators are not derivations 
and do not obey the chain rule for derivatives. We show how these problems may be resolved and illus- 
trate various uses of the characteristic. In particular, we establish the converse of Noether's Theorem 
for difference equations, we show (without taking a continuum limit) that the conservation laws in the 
infinite family generated by Rasin and Schiff are distinct, and we obtain all five-point conservation laws 
for the potential Lotka-Volterra equation. 



1. Introduction 

Current research in symmetry methods owes a tremendous debt to Peter Olver. In particular, his 
remarkable text, "Applications of Lie Groups to Differential Equations," remains pre-eminent after more 
than a quarter of a century. It is a masterpiece of scholarship that is notable for the lucidity of its 
exposition and the precision of its proofs. The first (1986) edition was the first text to describe the 
conditions under which the converse of Noether's Theorem |18j holds. A cornerstone of this result is 
the proof that, for any system of partial differential equations (PDEs) in Kovalevskaya form that is 
locally analytic, there is a bijection between equivalence classes of conservation laws and equivalence 
classes of characteristics. A simpler proof of this result, due to Alonso [m], is incorporated in the second 
edition [l9] of Giver's text. 

Here is a summary of the main definitions for scalar PDE^ For a given PDE, A = 0, a conservation 
law (GLaw) is a divergence expression that vanishes on solutions of the equation, so that 

DivF = when A = 0. 

A GLaw is trivial of the first kind if F vanishes on solutions of the PDE; it is trivial of the second kind if 
DivF = 0. A GLaw is trivial if it is a linear combination of the two kinds of trivial GLaws. Two GLaws 
are equivalent if and only if they differ by a trivial GLaw. If the PDE is totally nondegenerate (see [l9] ) 
- for instance, if it is in Kovalevskaya form - the GLaw can be integrated by parts to find an equivalent 
GLaw in characteristic form, that is, with 

DivF = QA. (1) 

The multiplier Q is called the characteristic of the GLaw. 
For example, the KdV equation, 

A = ut + + u^^x = 0, 

has a GLaw with 

DivF =Dt - ul^ + Qu'' + u'^u^^ ~ '^u^u^a^x + " 

= {u^ ~ 2u,D,) A. (2) 
Integration by parts yields the characteristic form of ^ : 

DivF = lu^ + 2u^^ A. 
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So Q = (v? + 2uxx) is the characteristic and ^ is equivalent to the CLaw 

A Qu^ - + A Qu'' + u^u^x + 2u^ut + ul^^ = 0. (3) 

A characteristic is said to be trivial if it vanishes on solutions of the PDE. By definition, the set of 
characteristics is a vector space; two characteristics are equivalent if they differ by a trivial characteristic. 
Therefore, the correspondence between equivalence classes of characteristics and CLaws makes it easy 
to identify when two seemingly different CLaws are equivalent: one only needs to compare their char- 
acteristics. Given a nontrivial characteristic, it is usually easy to reconstruct a corresponding CLaw by 
inspection; this can also be achieved systematically with the aid of a homotopy operator. In particular, 
where Noether's Theorem applies, each characteristic that arises from a one-parameter (local) Lie group 
of variational symmetries can be used to construct an associated CLaw. 

Until now, Alonso's result has not been transferred to difference equations. Yet one might wish to 
approximate a given PDE by a finite difference scheme that preserves difference analogues of several 
CLaws, particularly those that have a clear physical interpretatiorj^ This raises the question: is there a 
function that characterizes each equivalence class of difference CLaws? 

Of course, a given difference equation may be interesting in its own right, whether or not it is an 
approximation to a differential equation. Recent work has shown that CLaws of partial difference equa- 
tions (PAEs) have many features in common with CLaws of PDEs. For instance, Dorodnitsyn |5||6] has 
formulated a finite difference analogue of Noether's theorem. Hydon & Mansfield ^12j studied variational 
problems whose symmetry generators constitute an infinite-dimensional Lie algebra and derived a dif- 
ference analogue of Noether's Second Theorem. CLaws of a given PAE can be found directly (whether 
or not the PAE is an Euler-Lagrange equation); see 
PAE in Kovalevskaya form (see below) and see |22 
quad-graph equations. The main shortcoming of the direct construction method is that the algebraic 
complexity of the computations grows exponentially with the order of the CLaw; in practice, therefore, 
the method is restricted to low-order CLaws. 

Given a difference equation, A = 0, one cannot obtain the characteristics of its CLaw in the same way 
one does as for a PDE. For PDEs, the chain rule ensures that each CLaw is linear in the highest-order 
derivatives through which the dependence on A occurs. Integration by parts is then used to find the 
characteristic. The analogue of integration by parts for difference equations is summation by parts. 
However there is no analogue of the chain rule and therefore CLaws typically depend nonlinearly on A 
and its shifts. Consequently, it is not possible to construct a characteristic merely by summing by parts. 
In the current paper, we show how this difficulty can be surmounted. 

We also derive and use a difference analogue of Alonso's result. By pulling the characteristic back 
to a specified set of initial conditions, one can determine a function (the root) which labels the distinct 
equivalence classes of conservation laws. We show how the root is calculated in practice; examples 
include an integrable PAE with infinitely many CLaws. It is well-known that integrable PDEs have 
infinite hierarchies of CLaws which can be found using recursion operators, mastersymmetries, or Gardner 
transformations. Recently, Mikhailov and co-workers 16|17 and Rasin and co-workers f20'21' have shown 



To) for an algorithmic approach that works for any 
24 for applications of this approach to integrable 



that the same is true for integrable quad-graph equations. Having used the Gardner transformation to 



construct such a hierarchy, Rasin & Schiff 21 took a continuum limit in order to show that these CLaws 
are distinct. By using the difference analogue of Alonso's result, we show how to determine directly 
when CLaws are distinct, irrespective of whether they are preserved in any continuum limit or whether 
the underlying PAE is integrable. 

To establish the necessary results, it is helpful to begin by looking at scalar CAEs {^^- We define a 
characteristic of a first integral and show that the characteristic is trivial if and only if the first integral 
is trivial. In S|3) the definition of a characteristic is extended to CLaws of PAEs; we prove that there is 
a bijection between equivalence classes of CLaws and characteristics. This result has several immediate 
applications. In §4.1[ we use the characteristic to show that the CLaws in the infinite hierarchy for 
dpKdV generated by the Gardner transformation in |21 are distinct. Perhaps the most fundamental 
application is the establishment of the converse of Noether's Theorem (fj5]). Consequently, for Euler- 
Lagrange equations in Kovalevskaya form, there is a bijection between equivalence classes of variational 
symmetries and CLaws. Finally, we show how to use the characteristic to find CLaws of a given PAE 
(the potential Lotka-Volterra equation); this provides an alternative to the direct method. 



"^This is one of the oldest branches of geometric integration but, by exploiting the growing power of computer algebra 



systems, some new strategies for doing this have been developed recently jTj. 
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2. Scalar ordinary difference equations 

Although the main focus of this paper is on PAEs, it is instructive to look at scalar OAEs first. The 
independent variable is n € Z and the dependent variable is u € M. It is convenient to regard n as a 
free variable and to denote the shifts of u from a fixed but unspecified n by Ui := u{n + i). In order to 
evaluate first integrals on solutions of the OAE, one must be able to eliminate the highest (or lowest) 
shift of u. Therefore, we restrict attention to explicit K^^^-order OAEs, which are of the form 

A := - 7(n,Mo, . . . = 0, - — ^ 0, . (4) 

duo 

The set of initial conditions is the set of values z = {n, uq, . . . , uk-i} from which all Ui, i > K , can be 
calculated. 

A first integral of Q is a non-constant function, 0(n, uq, • . • , Wif-i), that is constant on solutions. It 
is helpful to introduce the forward shift operator, S'„, and the identity operator, /, which are defined by 

Sn ■■ {n, f{n), Uj) (n + 1, f{n + 1), u^+i), / : {n, f{n),Ui) i-^ {n, f{n), m); 

here / is any function that is defined at n and n + 1. In terms of these operators, (f> is constant on 
solutions if and only if the following difference CLaw holds: 

{Sn -l)(j) = when A = 0. (5) 

It is useful to refer to constant solutions of ^ as trivial first integrals, by analogy with trivial CLaws of 
the second kind. (Trivial CLaws of the first kind cannot occur when depends only on n, uq, . . . , uk-i-) 
A nontrivial first integral must depend on uk-i, otherwise (Sn — 1)4' does not depend on A (in which 
case, the only way for to be a first integral is to be identically constant). Therefore, the CLaw can be 
written as 

C(z, A) := 0(n+l,wi, . . . , uk-i, A + 7(71, wq, ■ • ■,uk-i)) - (l){n,uo, . . . ,uk-i), (6) 
where C(z, 0) = 0. We now use the Fundamental Theorem of Calculus to write the CLaw in the form 



C(z,A)=^^^^C(z,AA)dA 



/A=0 

=A / 0^K(n + l,ui, . . . AA + 7(n, uo, . . . rfA. 

(Throughout this paper, the partial derivative of a function, /, with respect to its i^^ continuous argument 
is denoted by f^i). By analogy with differential equations, we define the characteristic to be the multiplier 

{Sn<t))\uK=A+j - {Sn<l>)\uK=j 



dC{z,XA) ,^ 



As with differential equations, a trivial characteristic is one that vanishes on solutions, so that Q{z, 0) = 0. 

A trivial first integral is a constant, so C(z, A) = 0; therefore any trivial first integral has a trivial 
characteristic. To show that any trivial characteristic corresponds to a trivial first integral, it is helpful 
to define the root of the characteristic to be the function 



Q{z) :— lim (5(z,/iA) = lim 



('S'ri(?!>)|u;c=^tA+7 — {Sn4>)\uK=j 



M->o iiA 

^(f)^K{n + l,ui,..., UK-i,l{n, Uq,..., Uk-i)); (8) 

to do this, we require that Sn4> is differentiable in its K'^^ continuous argument at uk = 7- If the 
characteristic is trivial, its root is zero, so 

= (l>,Kin + l,ui, . . . ,UK-i,l{n,UQ, . . . ,uk-i))- (9) 

As 7^1 ^ 0, there is only one way for (j9]) to be satisfied identically in u^: the CLaw ^ cannot depend 
on A, so (/) must be a trivial first integral. 

Having dealt with the question of triviality, we now show how to reconstruct the first integral from 
the root. For clarity, we begin with a second-order example, but the same procedure applies in general. 
The OAE 

l:\:^U2-^{n,UQ,ui)=Q, 7(71, uq, ui) = - (^ui 5uo + 3 {l + (mi + -uq)^}^^^) , (10) 
has a first integral 

^(n, uo, ui) = 2-" [ui + uo + {1 + (mi + uo)^}^^^) • (H) 
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Consequently, the characteristic is 

Q(z, A) = 2-("+i) (a + {1 + (A + 7(n, uo, u,) + mf}'^^ - {l + {-f{n, Uq, ^i) + u^f}'^') /A, 
and so the root is 

g(z) = 2-("+i) (l + {7(n,uo,ui) + ui} {l + (7(n,uo,^^i)+ui)'}^'^') (12) 
22-" (5 + 2(wi + lio)^ + 2(ui + lio) {1 + (ui + uo)2}^/^) 



2' 

25 + 16(mi + uof 



(13) 



To reconstruct the first integral from the root, one must reverse this process. First use the OAE (10) 
to eliminate uq from (13), obtaining (12). Now treat ui and U2 as the continuous variables in (|8]), which 
amounts to 

^</>(n + 1, = 2-("+i) (1 + {U2 + «i) {1 + + ui)2}-'/') . 

Solving this and then applying gives 

cj>{n,uo,ui) = 2-" (ui + {1 + (mi +Mo)^}^^^) + /(?i,mo), 

where /(n, uq) is yet to be determined. The determining equation is ([5]), which amounts (after simplifi- 
cation) to 

f{n + 1, ui) - f(n, uo) = 2-("+i)(ui - 2uo). (14) 
Differentiating this with respect to Uq gives 

/,i(n,wo) = 2-", 

and so 

/(n,Uo) = 2-"uo +g(n). 



Thus, ((MJ) yields the OAE 

g{n + 1) — gin) = 0. 

Consequently, g(n) is an (irrelevant) arbitrary constant, which can be set to zero without loss of generality. 
This completes the reconstruction of the first integral (11). 

The process of reconstructing the first integral for a general scalar OAEs is similar. For any first 
integral, (j), of the i^*''-order OAE (|4|, the root is 

Q{z) := Q{z,0) = (l)^K{n + l,iti, . . . ,UK-i,-f{z)). (15) 

Given a root, eliminate uq in favour of uk to obtain 

— d 

Q{n,uo{n,ui, . . .,uk),ui, . . . ,uk^i) = 4>{n + 1,mi, . . .,uk-i,uk)- 

duK 

Integrating this and applying yields 

(j){n,uo, . . . ,UK-i) = / Q{n~ 1,uq, . . . ,UK-i)duK-i + f{n,UQ, . . . ,uk-2)- 



All that remains is to find /. Substitute (j) into ^ and simplify to obtain the determining equation 
for /. Differentiate this with respect to uq, then integrate to obtain / up to an arbitrary function 
g{n,ui, . . . ^uk-2)- Obtain the determining equation for apply S~^, and repeat the whole process 
with S~^g replacing /. Continue in the way until the remaining function to be determined depends on 
n only. The determining equation for this function (which we call h) is of the form 

h{n + l) - h{n) = H{n), 

where H{n) is given. The solution is obtained by summation; this completes the reconstruction of the 
CLaw. 
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Figure 1. A PAE in Kovalevskaya form: U30 = ^(to, n, Uq, Ui, U2). The box encloses 
the values Uij on which uj depends; these are represented by crosses. The dashed lines 
represent the initial conditions that would be required to obtain all Uij in the upper 
half-plane. 



3. Partial Difference Equations 

We now generalize the ideas from the last section to scalar PAEs for u G M with two independent 
variable^ n = {m,n) £ 1? . Again, we regard the independent variables as being free; given {m,n), let 



u 



u{m + + j). [The indices i and j may be negative: each minus sign in the subscript should 
be treated as being attached to the following digit. For instance, ui_i denotes u{m + l,n — 1).] The 
action of the shift and identity operators on {m,n) induces an action on every function f{m,n), and in 
particular on u^j, as follows: 



I 



{m,n,f{m,n),Uij) ^ (to + 1, n, /(m + 1, n), 

(to, 71, /(to, n), U,j) 1-^ (to, 71 + 1, /(to, 71 + 1), + 

(to, 71, /(to, 71), ) (to, 71, /(m, n),Uij). 



A PAE is written as 



A(to,7i, [m]) = 0, (16) 

where [•] denotes the argument • and a finite number of its shifts. 

A CLaw of a PAE is a divergence expressior{^ that vanishes on solutions of the PDF: 

DivF := {Sm ~ I)F + {Sn - /)G = when [A] = 0. (17) 

The functions F := F{ni, n, [u]) and G :— G{m, n, [u]) are the densities of the CLaw. In the same way 
as for PDEs, a CLaw of a PAE is trivial if and only if it is a linear combination of the following two 
kinds of trivial CLaws. 

First kind: F|[a]=o = 0. The densities vanish on solutions, so we call these trivial densities. 
Second kind: DivF = 0, without reference to the equation A = and its shifts. This occurs if 
there exists a function H such that F — {Sn — I)H and G = —{S„i — I)H . 

3.1. Kovalevskaya Form. A crucial step in dealing with first integrals is to replace uk by A + 7. We 
can do something similar for a /^th-order PAE ( 16 1 if it is in Kovalevskaya fornrA 



A := UA'o - '^("i, 71, uo, ui, . . . , uk-i) = 0, (18) 

where Uj = {uij : j G Z} and there exists j such that duj/duQj ^ 0. A schematic example of a 2-D scalar 
PAE in Kovalevskaya form is shown in Figure [T] Let z = {m, ri, Uq, Ui, . . . , Uk-i} be the (minimal) 



initial conditions from which shifts of ( 18 ) can be used to find any point of the form uij for I > K . The 
function uj := uj{z) depends on a finite subset of these points. It is convenient to denote shifts of w by 
ujij := S'^Siiv. 

A scalar PAE with two independent variables is explicit if it can be transformed into Kovalevskaya 
form by an admissible change of independent variables, that is, by a bijective linear map from to 
itself. The new independent variables are 

^j^A(^j+h where A e GL{2, Z), det(^) = ±1, and b e Z^. (19) 



''The corresponding results for systems of difference equations with arbitrarily many independent variables are obtained 
mutatis mutandis; see jrj. 

■^For differential equati ons on and difference equations on Z^, the set of divergence expressions is the kernel of the 
Euler-Lagrange operator [l3| . 

'^For a given PDE in Kovalevskaya form, any equation that holds on solutions of the PDF can be pulled back to an 
identity on the initial conditions; with the above definition, the same is true for PAEs. 



6 



TIMOTHY J. GRANT AND PETER E. HYDON 



U21 




Figure 2. Transformation of a quad-graph equation into Kovalevskaya form. The 
dashed hnes show the initial conditions that would be required to determine all Uij 
(respectively iiij) in the upper- right (respectively upper) half-plane. 



Although the value of the u at each point is unchanged, the coordinates of the point have changed, so it 
is helpful to define 

u{m,h) :^ u{ni{m,n),n{rh,h)). (20) 
Now fix (m, n); using the shorthand — u{m + i, n + j) and setting uqo = itoo, we obtain 

Uij = Sl^Siuoo = u{m{rn + i,n + j),h{m + i,n + j)) = U{m{rn+i,n+j)-7nini,n)}{n{Tn+i,n+j)-n{rn,n)}- 

For instance, the shear 

A^a ]). b.(») (21) 



transforms any quad-graph equation, 

uii = w(to, n, uoo, uio, uoi), 

into the Kovalevskaya form 

U21 = U!{m{m, n), n{rh, n), uoo, Wio, uii). 

In particular, the dpKdV equation (HI in the ABS classification [l]), 

P — a , , 

uu = uqo H , (22) 

uio - uoi 

is transformed into 

(3 — a , , 

"21 = -"oo + — ■ (23) 



Lemma 3.1. When an explicit scalar PAE is transformed according to ( 19 ) and [20], there is a bijective 



correspondence between equivalence classes of CLaws of the original PAE and the transformed PAE. 
Proof. The transformation is of the form 

^) = (c + ^d-bc^±l, a,b,c,d,eJeZ, 

so the effect of the original shift operators on the transformed independent variables is 
S,nrh — a(jn + 1) + 6n + e = m + a, Smfi = c(m + \) + dn -\- f — fi + 
Snrh ~ am + b{n + 1) + e = m + &, Snfi = cm + d{n + 1) + f = n + d. 

Thus 

Q na nc n nb od 

Therefore given a CLaw with densities F and G, 

{S^ - I)F + {Sn - I)G - iS%,Sl - I)F + {Sl,Si - I)G 

= [{S^^ - I)F + {Si - I)g] + [{SI - I)S^^,F + {Si - I)SiG} 

= {S,n - I)F + {Sn - I)G, (24) 

where F and G denote the densities in terms of the transformed variables and the last equality is 
obtained by factorizing each expression in braces. If the original CLaw is trivial of the second kind then 
it vanishes identically, so the transformed CLaw must also vanish identically. If the original densities F 
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A = (-l)'"{2«oo{«ll -«lo) + a-^) 
G, = i-ir (2«io (no-1 + .n^^) - 

-P2 = («00 - «10) (m00«10 - a) - (UQO - Mil) {«oo"ii - /3) 
= ("10 - no-1 - ,,^^) («10 («0-i + ,,^^) ^ 

F3 = (-l)™ («oo("ll - "10)(«00 + Mio + Mil) + a(MOO + Uio) - /3(moo + Mil)) 
G3 = (-1)- (mio + Mo-1 + .r^^) (mio (mo-1 + ^^^) - a) 

Fi = (-l)'"+i (2moo^(«io^ - + 4moo(/3mii - auio) + - /jZ) 

64 = (-1)'" (2mi„2 (mo-1 + ^^^^f^Y - 4aMio (mo-i + ^.^^f^) + 

^5 = - 1- (=lf5^) + 1- ("0-1 - "00 + ,,^^) 

= In (mi_i - Mio + (/3 ~ a) («o-2 - «o-i + - ^7^^) 

Fg = - In f Moo - "0-1 + ) + In f ) 

Ge = In f (/3 - a) f mo_2 - «0-i + ^^^^^^ fciS^) "^^ 

Ft = [m- n)F3 + nFe 

(57 = (m - n + 1)65 + nGe 



and G are trivial then so are F and G; consequently, the new densities F and G will also be trivial. As 
the transformation is invertible, the converse is also true, which completes the proof. □ 



For example, when a quad-graph equation is transformed by (21 ), the transformed equation has CLaws 
with the densities 

F = F + G, G = S,nG. (25) 
To transform back to the quad-graph, use 



m \ I 1 —1 \ j m 
n J [ 1 [ h 



(26) 



Then the densities for the CLaws are 

F — F — Srf^G^ G — Srf^G. (2'^) 

Generally speaking, the new densities will depend on variables other than the transformed initial condi- 
tions; however, the difference equation can be used to pull them back to expressions that depend only 
on the transformed initial conditions. 

For later use, the dpKdV equation has four three-point and three five-point CLaws [24]; the trans- 
formed equation ( 23 1 has corresponding CLaws whose densities are listed in Table [l] (where the tildes 
have been dropped to prevent clutter). 

3.2. The characteristic. The characteristic of a given CLaw of a PAE in Kovalevskaya form (18) 
is defined as follows. Any Uij may be written in terms of the initial conditions, z, and [A]. Given 
a differentiable function, G(z, [A]), that satisfies G(z, [0]) — 0, the Fundamental Theorem of Calculus 
yields 

G(z,[A])=/' -^G(z,[AA])dA= /' V (^^^^^A)^^^^^^!^ dA. 
Summing by parts gives 

G(z, [A]) -A {E^ (G(z, [A]))}|^^^^ d\ + DivF, (28) 
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where _Ba is the difference Euler-Lagrange operator that corresponds to variations in A, 



E^{C{^, [A])) := y S-:S-^ , (29) 

and where F has two components, each of which vanishes on solutions of the PAE. (The generahzation 
to PAEs with more than two independent variables is obvious.) If, in addition, C(z, [A]) is a divergence 
expression (so that it is a CLaw) then, by subtracting the trivial CLaw DivF from both sides of (281, 
one obtains the equivalent CLaw 

C(z, [A]) - C(z, [A]) - DivF = Q(z, [A]) • A 

where 

Q(z,[A]):=/ {(i?A(C(z,[A])))}|^^^^dA. (30) 
J\=o 

The function Q is a multiplier of the difference equation, so to be consistent with the continuous case 
and 15 , it is defined to be the characteristic of the CLaw. 

Just as for OAEs, a trivial characteristic is one that vanishes on solutions, so Q(z, [0]) = 0. To factor 
out trivial characteristics, we define the root of Q to be 

Q(z) hm Q(z, [^.A]) - hm {i?^A(C(z, [p^A]))} \^^^^dX = {E^{C{z, [A]))} \^^^^^. (31) 

For PDEs, the characteristic of a CLaw is trivial if and only if the CLaw is trivial. We will now show 
that, with our definition of the characteristic (30), the same is true for PAEs in Kovalevskaya form. 
Hence, we need to prove that a CLaw is trivial if and only if Q{z) ~ 0. 

3.3. A trivial CLavif implies a trivial characteristic. A trivial CLaw of the second kind (DivF = 0) 
vanishes identically, so the proof is immediate. However the first kind of triviality takes the form 
F(z, [0]) = 0. To deal with this wc use the identity 



Ea (DivF(z, [A])) - EA{iS„, - I){F{z, [A]) - F(z, [0])) + (S,, - /)(G(z, [A]) - G(z, [0]))^ 

' dF{z,[A]) ^ dG{z,[A]) 
dSl^SiAi d{Sl^SiA) 



dSmZk ^ d{F{z, [A])~F(z,[0])}^ _ ./ 



I ^^-,^-J^ dSrnSl^SiiA dF{z,[A]) dSnSLS^^A dG{z, [A]) \ 

where e z and we have used the fact that (for equations in Kovalevskaya form) SnZk G z. The identity 
. , ( dSmSl,Si;A dFjz, [A]) \ >^ 1 ; dF{z,[A]) 



together with a similar identity for the G term, simplifies ( 32 ) to 

£;A(DivF(z, [A])) = J2 ^^^n' ( ffZ7^ ■ [^]) - M)) 

From this, it immediately follows that 

{i?A(DivF(z,[A]))} 

so the root is zero. 



[A]=0 



3.4. A trivial characteristic implies a trivial CLaw. Having proved that the characteristic of a 
trivial CLaw is a trivial characteristic, it is now clear that trivial densities may be added to the densities 
of any CLaw without affecting the root. In particular, the components of F(z, [A]) — F(z, [0]) are trivial 
densities, so any explicit dependence on [A] in the densities can be removed. Therefore, we now assume 
that the densities of any particular CLaw for ( 18 ) do not depend explicitly on [A] ; with this assumption, 
only the second kind of triviality can occur. Consequently, the CLaw has densities of the form 

F := F(m, n, Uo, . . . , uk-i), and G := G(m, n, Uq, . . . , uk-i). 



DIFFERENCE CONSERVATION LAW CHARACTERISTICS 



9 



S,„F 



*■ U2 
►-Ui 



Figure 3. A graphical representation of the terms in the characteristic for a scalar 
PAE with two independent variables. The solid black box encloses points on which 
SmF depends and the dash-dot box encloses points on which S~^S„iF depends. The 
dashed boxes enclose points on which luo-r, wqo and (jJqr depend. 



Clearly, SnG does not depend on Uk, so the root is 
_ ^ . g 

Here we have assumed that F depends on {uK-i.j,j = 0, . . . ,R} and on no other uk^ij', there always 
exists an i? > for which this assumption is valid, as we have the freedom to select the starting-point 
(TO,n) relative to which all Uij are compared. The shifts of SmF on which Q depends are shown in 
Figure [sj here and henceforth (for brevity) , we refer to the values Uij that occur in any expression as 
'points' on which the expression depends. 

Let mql be the leftmost point in Uq on which uj — wqo depends. Then 

- — °^^0, jeZ; - — ^ = 0, i>j ijez. 

C"0(L+j) CUo(L+j) 

On solutions of the PAE, each point Ugii^jj (such as the square points in Figure [s]) may be replaced by 
utQj (represented by the discs in Figure h5| as an independent variable in Q, because the determinant of 
the Jacobian 5("o-fi.-v"oH) — ^ jg nonzero. Therefore 

R Q 

Q = y] ^n'^ F{m + 1, n, ui, . . . , UK-i, KD EUSrr.F). (33) 

This is a restricted difference Euler operator corresponding to variations in w; in effect, it treats m and 
the Ui terms as parameters. The kernel of this operator is made up of total divergences of the form 
{Sn — I)H, and functions of z\uo only (see Lemma 3.2 below). Thus if the characteristic is zero on 
solutions of the PAE then 

F(TO-|-l,n, ui, . . . ,uk-i, [ujqj]) = (Sn - I)H{m+l,n,Ui, . . . ,uk-i, [ujoj]) + f{in+l,n,Ui,. . . ,uk-i), 
for some /, and so 

F(z) = {Sn - I)H{m, n, uo, . . . , uk-2, [u(^k-i)j]) + fim, n, Uq, . . . , uk-2)- 
Adding the trivial CLaw 

FT^-{Sn-I)H, GT^{Sm~I)H, (34) 

to the original densities gives the equivalent densities 

F = /(m, n, uo, . . . , UK-2), G = G+ (5,„ - I)H. 

The G density may contain A terms, but these can be removed by adding a trivial density. Thus the 
divergence expression for these densities cannot contain any A terms, so in order for it to be a CLaw it 
must vanish identically and is thus trivial. 

The following lemma identifies the kernel of the restricted Euler operator; part of the proof will be 
used in §6 to enable us to reconstruct a CLaw from its root. 

Lemma 3.2. The kernel of E^^ consists of sums of functions that are independent of uj and its shifts, 
together with total divergences in the n direction. 

Proof. In the following, we use the notation u — {uij : 1 < i < K — 1} and 

uj\ := XujQo + (1 - X)gim, n, u). 
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Table 2. Roots of the transformed dpKdV 



n 




Q2 


= (itii — uio)(ttii + uio — 2a;oo) + a — /3 


Qz 


= (—1)™" ((uio — + Mio + 2tj()o) + a — /3) 


O4 


= 4 (-1)™ (0^00(^^10^ - Itii^) + Qliii - /? -Uio) 




= (/? a)(a;o-i cjoo) (itio + u^o-i-^ot^ q) (^^oo ojoi) (^^^ii ?ti2 + J 
+ (cjoo — ^01)""^ — (1^0-1 — i^ooY^ 




= {13 a){ojoo ujoi) (uii 'iio + J (/3 a) (a^o-i woo) (uio + 

+ {u!o-l - OJoo)"^ - (i^oo - ^Ol)"^ 



Q7 = W~ ")(woo - a;oi) ' I (uii - ui2 + ) + (^11 ^ "10 + J I 



+ {rn — n) Qg + nQg + 2 (cjoi — ^oo) 



where g is any convenient function (usually, g = 0). For any differentiable function / = f{m, n, u, [uj]), 
-^/(TO,n,u, [ujx]) = ^ni^oo - g(n, u)) 

= (woo-.9)-£^cja(/('^)'^iU7 [i^a])) + (.Sn - I)h{m,n,u, [w], A), (35) 



for some function h. Integrating (35) with respect to A, we obtain 

f{m,n,u,[uj]) = f{m,n,u,[g]) + {uJoo-g) / E^^{f{m,n,u,[ujx])) d\+{S„~I) / h{m,n,u,[uj], X) dX 
If / e ker(iJ;^) then the result follows. □ 

4. Using roots to detect equivalence 

For PAEs, the fact that a CLaw may involve a large number of points can make it difficult to 
identify its underlying order. The results of §3 have established that the root completely characterizes 
an equivalence class of conservation laws, resolving this difficulty. 

As an illustration, consider the transformed dpKdV equation (23). Using maple, we have calculated 
the roots of the CLaws of ([23| that are listed in Table [I] These roots are displayed in Table [2|^ expressed 
in terms of the functions (because this is more compact than pulling back to write each in terms 
of z). The transformed dpKdV equation also has the following CLaw 

F = (-1)"^+^ (2 uii ( uoi + ) -13]+ i-ir^^ (a - 2un ( uoo + ^ ^ 





— a 


f 11 


- U12 




— a 




— till 




— a 



UlO — I'll 



G = Moo H Uoo H 1*0-1 H + 



y \\ Uio — Itll / \ Mio— Ull/ y U\-\ — U\o , 

(-1)"+^ / / P-a \ ( P-'^ \ oaf 



ce Uoo H + " "0-1 ^ 2/3 Uoo + 



where y =uoo — uo-i 



y V \ "10 ^ "11 / \ "1-1 ^ "10 / \ ^10 — uii 

13 -a (3 -a 



Ui_i — Uio Itio — Uii 



This apparently high-order CLaw's root is 2 (—1)™ (un — uio), which shows that actually it is equivalent 
to a multiple of the first CLaw in Table [T] 
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Figure 4. The figure on the left shows the densities for the third CLaw in the hierarchy; 
on the right the extreme shifts of the densities in the characteristic are shown. 



4.1. The Gardner CLaws for dpKdV. In 21 , Rasin and Schiflf used a discrete version of the Gardner 
transformation to construct an infinite number of CLaws for the dpKdV equation (22). (Rasin f^O" has 



subsequently used the same approach to generate an infinite hierarchy of CLaws for all the equations in 
the ABS classification and one asymmetric equation.) They showed that these CLaws were distinct by 
taking a continuum limit and showing that the resulting CLaws for the continuous equation are distinct. 
By using roots, one can prove that their CLaws are distinct without having to take a continuum limit. 

In this section, we consider the dpKdV equation on the quad-graph rather than in Kovalevskaya forn|^ 
The dpKdV equation can be solved for any of the points on the quad-graph; hence we choose the initial 
conditions z = {m, n, Uio, u_ji, Mofe|«, J, G N, A; € Z} which are shown by dashed lines in Figure [4] 

The densities for the CLaws generated by the Gardner transformation are the functions Fi and Gi in 
the expansion of 

oo \ oo 

1 + V 4o = Y] F^^\ (36) 

/ ^ OO \ oo 

G=lne-lnKo-«2o)+ln l + ^V4'''^eM =lne + VG,e\ (37) 



V, 



00 i=l / i=0 



in powers of e; here 



^00 — „ „ ' "oo — „, „ 00 "lO ' 

"00 — "20 "00 ~ "20 , 

J = l 

and iiqq is referred to as the i*"^ order v term. To prove that the CLaws are distinct the following lemma 
is used. 

Lemma 4.1. For each a E N, 

^-^4^V0 and ^^4^)=0, j>2. (38) 

C'"(o!-|-l)0 O'"(a+i)0 

Proof. Proof is by induction; the base case {a — 1) is immediate. Assume that the lemma holds for 
a = k — 1, where k > 2. For each i £ N, 



^ (fc) 1 I \^ f U}\ (k-f) , U) 9 
'^oo = a ^00 no 



The first term in the summation is zero for all j, because the highest value j can take is /c — 1 but the 
lowest value of i is 1. Similarly, the second term vanishes for all j > 2, so only one term is left, namely 



d (fe) / (ii\2 ^ f d 
du ""00 



Using the induction hypothesis, if i = 1 then (39 1 is nonzero and if i > 1, (39) is zero. □ 



A consequence of this lemma is that, for a > 1 



"^^^^'^ ^-"^^"la^I^^oo j^O. (40) 

^Kovelevskaya form is convenient for proving that the root characterizes each equivalence class of CLaws. Hovfever, 
for any explicit PAE, the root (with respect to an appropriate set of initial conditions) can also be calculated without 
transforming to Kovalevskaya form. 
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As the second factor in (40) doesn't depend on uqq, we obtain 

„(" + !) \ 



"00 



a+2)0 



„(1) 



^0. 



Expanding out (37l shows that, for a > 1, the highest order v term in Ga is ^^oo^^V^oo''' 



(41) 



by Lemma 



this is the only term in Ga to depend on U(q+2)o- Therefore S'„(wqq /"yoo ) is the only term in 
the CLaw to depend on M(q,+2)i and so to depend on S'J^+^A (because SmF depends only on points in z 
except for un - see Figure III). Thus the root is 



Ei^^SmFa + ^nGa)\!^^_Q — 



dSmFa 

dA 



E ^ dSnGg c<-(«+l)c f 9 
1=0 '^^rn^ \OU(a+2)0 



(q+1)- 
-00 

''oo 



[A]=0 



The final term is the only one to depend on u_^a+i),i\ no other term is shifted as far back. From (41 1, 



-(a+l)l 



-E^{S„,Fa + 5„G„)|rAi=o = S-^'^+^^S.. 



duoodU[a+2)0 



^o 

"oo , 



^0, 



so this term cannot be a linear combination of the other terms. Therefore, the characteristic does not 
vanish on solutions of dpKdV and so the CLaw is nontrivial. The roots of the lower-order CLaws with 
densities (i^i, Gi), . . . , (i^Q_i, Gq,_i) do not depend on u_(o,+i)i, so the CLaw {Fa, Ga) cannot be a linear 
combination of these CLaws and their shifts. Thus the CLaws generated by the Gardner transformation 
are distinct. 



5. The converse of Noether's Theorem 

One of the most useful significant applications of the root is to establish that the converse of Noether's 
Theorem holds for difference equations. Given a Lagrangian, L, the Euler-Lagrange equation is 

Each symmetry generator for the Euler-Lagrange equation is the prolongation of 

X = q/- ; (42) 

the function Q is the characteristic of the symmetry generator. In particular, X generates variational 
symmetries if XL is a total divergence, in which case 

E{XL) = 0. 

In this case, there exist functions / and g such that 

(SlSiQ) - (5™ - /)/ + (5„ - I)g. (43) 



Summation by parts is then used to rewrite ( 43 ) as 



Q ■ E{L) = E QS^^S-'^ l^-^L ) = {Sm - I)F + (5„ - /)G, (44) 



for functions F and G whose precise form is irrelevant. Thus if Q is the characteristic of a variational 
symmetry generator, it is also the characteristic of a CLaw for the Euler-Lagrange equation. Two 
variational symmetries are equivalent if they differ by a symmetry whose characteristic vanishes on 
solutions of the Euler-Lagrange equation (i.e. a trivial symmetry). Therefore, if the Euler-Lagrange 
equation is explicit (in which case the CLaw is trivial if and only only if the root is zero), there is a 
bijective correspondence between equivalence classes of variational symmetries and CLaws. 
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6. Reconstruction of CLaws from roots 

If the characteristic of a CLaw is known then the densities for the CLaw can, in principle, be re- 
constructed using homotopy operators [IT 19 . For PDEs (given an initiahzation) , the root of a CLaw 
can be calculated by the same approach as we have used; the root is necessarily a characteristic as a 
consequence of the chain rule. By contrast, the root of a CLaw for a PAE may not be a characteristic. 
So the key step in reconstructing a CLaw from its root is to find a characteristic which has that root. 



Our starting-point is the proof of Lemma 3.2 Replacing / by Sm.F and using the definition (33), we 
obtain 



F{m + 1, n, u, [w]) = (wqo - g) Q{m, n, u, [ujx]) dX + (5„ - I)H{m, n, u, 

Jx=o 



M) + /(?7i+ l,n,u), 



for some / to be determined. The H term can be set to zero without loss of generality by adding a 
trivial CLaw of the second kind, so we can assume that 



F(m-M,n, u, [w]) = (wgo-g") / (3(m, n, u, [wa]) rfA -I- /(m -I- 1, n, u). 



(45) 



In general, Q, and as a result (45), contains negative shifts of uj. These can be removed term-by-term 



by adding trivial CLaws of the second kind, until one obtains an equivalent density SmF that has no 
negative shifts of lo. 

Equation ([45| (shifted if necessary, as discussed above) contains all of the w-dependence of the CLaw. 



Having obtained this, a characteristic for the CLaw is calculated by replacing ujqj in (45) by 5*;^ A + ujQj. 
From dSOl), 



g(z, [A])= I i;A(F(m-fl,n,u, [A-f a;oo]))|AH>AAdA, 



(46) 



and the CLaw can be written as 



C := Q{z, [A]) • A Q(z, [uko - ^^(z)]) {uko ~ c^(z)) • 

Homotopy operators (see [IT]) may then be used to find the densities. Alternatively, once the w depen- 
dence has been found, the arbitrary function / and the other densities can be constructed directly by a 



variant of the method that we used for OAEs (see 10 for details). 



6.1. Example: reconstruction of CLaws of dpKdV. To illustrate this, we will reconstruct two 
CLaws of dpKdV from their roots. For both examples, we choose g{n,u) = 0, so uj\ = Xojqo- (In 
practice, we start with this choice and only change it if the integral is singular.) First, we use the fourth 
root in Table [2] which gives 

F{m + 1, n,uio, uii,u}oo) = cjoo / 4 (-l)'"+^{(uio)^Aa;oo - XuJooiuu)'^ + aun- Pmo} dX + f{m + l, n, uio, un) 
= 2 (— l)™+"^a;oo {t^oo ((itio)^ — (uiif) +2(auii — Pmo)} + f {m+ l,n,uio, un) . 

(47) 



This has no negative shifts of uj so, from (|45|) and (|46|), a characteristic is 

Q(Z,A): 



d 



A^O^A 



{2 (-1)"+^(A + ojoo) ((A + c^oo)("io' - un^)+2{aun-Puio))} 



Ah^AA 



=2 (-1)™+^((A + 2c^oo)(wio' - uii') + 2(aMii-^uio)) 



Thus 



C -.^QA = 2(-l)™+^ (m2i - woo)((m21 +woo)(uio^ - uii^) + 2(q!Uii-/3uio)) 

=2 {-l)"'(l5 + a- U21U10-U00M10-M21M11-M00M11) iu2iUio-U2iUii-uooUio+uooUii-l3+a) . (48) 
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This is a total divergence, so applying the homotopy operator from gives the required densities 

1 

6 



F = i (-1)™+^ (4uioMo-i/3 + 20-UOO/3 "ii + 4u_iiauoi + 2u_ioUoi/3 - lOuoo^wn^) 



+ ^ (-1)"^^ (uoo^u_i_i^ - 2uoo/3u-i-i - 4uooU-ioa - lOuooUioa + 2u^io^-uoo^) 

+ ^ (-1)"^^ (7moi^"ii^ - 14aMiiWoi - u-io^uoi^ - 2uio^uo-i^ - 2uoi^w-ii^ + Suqo^wio^) , 

G = ^ (-1)"' (-IOW10M20Q: + 2m_io^uoo^ - wio^uo-i^ - 5uoo^?iio^) 

+ ^ (-1)" (-2uoo/3w_i_i + lOwoo^^ioa + uqo^U-i-i^ + 2uioUo-i/3) 

+ TT (4w2oWi-i/3 + 5uio^U20^-2u2o^Mi-i^-4uooU-ioa-12n/3^ + 12na^) . 

6 ^ ' 

These densities contain points of the form w-ij and U2j- To find equivalent densities that are given solely 
in terms of z, one must shift the CLaw forwards and then use the dpKdV equation ( 23 ) to pull all terms 
back onto the initial conditions; this leads to even longer expressions! 

In practice it is much easier to use the direct construction method, which leads to more compact 
expressions for the densities. Starting from (47), we know that the densities are of the form 

F =2 (-1)'" uii(uii(uoo^ - uoi^)+ 2 (auoi-/3woo))+/ (m, n, uqo, uoi), 
G ^G{m,n,uaQ,uiQ). 

Substituting these into the CLaw and evaluating the result on solutions gives 



= C|a=o =2(-1)"^^ {uoo + — {uoo + — Ko'-wii')+2(aun-/3uio) 

- 2 (-1)™mii(mii(woo^ - uqi^)+ 2 (auoi -/Suqo)) + / (m + mio,uii) 

- / (to, n, uoo, woi) + G (m, n + 1, mqi, un) - G (to, n, uqo, uio) ■ 
Differentiating this expression, we obtain 







duooduoi duooduoi 



Consequently, / — f{m,n,uoo) + /i(m,n,uoi); however, h can be set to zero by adding a trivial CLaw. 
Therefore 

0=— =4(-l) MiiMoi+4(-l) uiia + - G (to,7i + l,Moi,uii) , 

duoi auQi 

and so 

G (to, n. Moo, uio) = 2 (-1)'"+^ uio^uoo^ + 4 (-1)™ uioa uqo + 9 (m, n, uio) . 
Dropping the tilde, 

Q ^ ^GIa^o ^ df 
which yields / — f{m, n). The final differentiation is 

so 5 = g{m,n). The CLaw now simplifies to 

{S„, - I)f{m, n) + (5„ - /).9(m, 7^) = 2 (-1)" (a^ - /S^), 
a solution of which is / = and g ~ 2n(— 1)™(q!^ — /3^). So the reconstructed densities are 

F =2uii (-1)" (uii(moo^ - "01^) + 2aum - 2/3 uqo) , 

G =2 (-1)" (-uio'uoo' + 2Mioauoo + '^(a' - /3')) , 
which are equivalent to the densities found by the homotopy method, as they have the same root. 
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For a more complicated example, consider the densities of the sixth Claw in Table [T] Its root depends 
on cjo-i, Woo and wqi- The characteristic is constructed so that terms depending on wq-i do not depend 
on uiQi- The term that depends on wq-i is 

-1 



/ii(m,n,ui, wo-ijWoo) := (wo-i-woo) - (/3-a)(i:Jo-i -t^oo) "lo H ui-i 

Then let 

/i2(m, n, Ui,woo,woi) ■■=ojQoQe{m,n,Ui, [\uj]) + (5„ - /) (woo^iiCn, Ui, [Aw])) 

^00^10 - t^OlUlO - a^OQMll + WpiMll 

-Awoo^^ii + Awoouio - AwoiMio - /3 + AwoiUn + a ' 
which is a function that only depends on woo and cjoi; this is integrated with respect to A to obtain 



'S'm^|[A]=o = y" /i2 dA + /(m + 1, n, uiq, un, M12) 



In I ^ _ ^ ) + f[rn + 1, n, itio, Uii, iti2)- (49) 



As in the last example, we could calculate the characteristic from (491 and use the homotopy operator 
to find densities for the CLaw. Once again, however, the direct construction method is preferable. By 
shifting (49) and choosing G to depend on appropriate values of z, the densities have the form 

, f {ui2 - uii){uoi - uqo) + a - f3\ 
^ = -ln| ^ _ ^ 1+/ (m, n,uoo,woi,uo2), 

G ^G {m,n,uoo,Uio,uoi,uii) . 



The dependence on the dpKdV equation is already determined by ( 49 ) , so all that remains is to find / 
and G. The same process (differential elimination followed by integration) is used as before. Skipping 
the details, we obtain 

, /(ui2-uii)(uoi-uoo)+a-/3\ w ^ 

F=-lnl ^ _ ^ I, G = -ln(uio - uii) , 

as required. 

7. Finding CLaws 

7.1. The Adjoint of the Linearized Symmetry Operator. The Gateaux derivative of a functional 
P is the operator defined in [IT] by 



^ .^s , , P \u + eQ\u\] - P\u\\ {d , 
Dp(Q) ^ lim ^ \^\—P[u + eQ[u]] 



de 

Explicitly, the Gateaux derivative of P is the shift operator with entries 

ij 

Therefore its adjoint with respect to the ^2 inner product is 

du 



e=0 



Ti* — { g-' Q-.i 1 Q-iq-j 



and so the Euler operator, E, is defined by 

E{P[u])=B*p{l), 

just as for PDEs. Using the Leibniz rule, 

EiP-Q)^D*p.Q{l)^D*p{Q) + D*QiP). (50) 
The action of the vector field X = Qd/duQQ on the functional P is 

pr XiP) = Dp(Q). 

So the linearized symmetry condition for a given difference equation, A = 0, is 

0= Da(Q)|[^]=o. 
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Table 3. Solutions of the ALSC for the potential Lotka-Volterra equation 



Q2 = 



Q 



M0-1«00 Moo M00^("10+"0l) 



I f02t'10 "10 

"00 moo^(mio+"oi) Uo-lUoO 

m,MioM-io (m+l)M20 

"oi"oo("io+«-io) «oo(«20+"oo)(mio+moi) 



4 UoiUooC'llO+U-lo) tiOo(M20+«Oo)(MlO+MOl) 

_ (-1)'" + "(M10^«01+U00^(f 10+«0l)) 

^5 tioo^"oi(«io+"oi) 



"10 "01 — "00 ("lO+MQl) 
"00^"0l("10+tiOl) 



The Euler operator acting on a expression is zero if and only if that expression is a total divergence [llj . 
Therefore Q is a characteristic of a CLaw if and only if 

o = £;(g.A) = Dl(g) + D^(A). 

Restricting this to solutions of the difference equation gives a necessary condition for Q to be a charac- 
teristic: 

0^-D*Am[A]=o- (51) 
In other words, the characteristics are members of the kernel of the adjoint of the linearized symmetry 
condition (ALSC), restricted to solutions. Arriola ^ showed that (51) must be satisfied for first integrals 



of autonomous ordinary difference equations. In a notable paper that introduces the idea of co-recursion 



operators for integrable difference equations 17 , Mikhailov et al. define a cosymmetry of a difference 
equation as being a member of the kernel of the ALSC. They state (parenthetically) that cosymmetries 
are characteristics of CLaws, but do not justify this. 

In this section, we give an example of a PAE where not every cosymmetry is a characteristic of a 
CLaw (see [2j[3] for a discussion of this point for differential equations). First, we find functions Q that 
satisfy ( |51[ ), using methods similar to those used to find symmetries of difference equations. Then we find 
additional constraints on Q by applying the difference Euler operator to QA. Finally, we reconstruct the 
densities using homotopy operators or inspection. For brevity, we omit most details of the calculations. 

7.2. CLaws of the Potential Lotka-Volterra equation. The potential Lotka-Volterra equation 
(pLV) 

— — 1. (52) 

uoo Uio 

is an integrable equation on the quad-graph. (It belongs to Class 4 of Hietarinta & Viallet's classification 
of quadratic quad-graph equations with polynomial degree growth 8|, and is a potential form of the 
discrete Lotka-Volterra equation introduced by Hirota & Tsujimoto |9 .) 



Rasin & Hydon's method for finding symmetries of quad-graph equations 25 is readily adapted to 
find solutions of the ALSC, shifted (for convenience) to 

= {{Q - Sn (w,2(5) - Sra (t^.sQ) " SmSn {u!,iQ))} |[^]^q- (53) 

We search for solutions of ( |53[ ) which are pulled back onto the initial conditions z ~ {m,n,UiQ,UQj}. In 
particular, we will look for roots of the form 

Q = Q{m, 71, u_io, wo-i, Uoo, uio, ""oi, U20, ^02)- (54) 

The corresponding CLaws depend on m, n, 7i_io, uo_i, Mqo, and Uqi only, and are therefore called 
'five-point CLaws'. 

By differential elimination and integration, one obtains the solutions of the ALSC; these are listed 
in Table [3j However, when one tries to use these roots to construct the corresponding CLaws of the 
pLV equation, the algorithm fails for (see Table It turns out that E{QA) ^ 0, and therefore 



this cosymmetry does not correspond to the characteristic of a CLaw. This is not surprising, as (51) is 
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The solution of the ALSC 

p _ (■"-io-'"oi)("o-i + "-io) 


is not the root of a characteristic. 
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necessary but not sufficient; nevertheless, to the best of our knowledge, this is the first example of such 
a cosymmetry of a PAE. It raises an interesting question: can an additional constraint be found that 
guarantees a solution of the ALSC is a root of a CLaw without the need to work through the (lengthy) 



process of reconstructing the characteristic? For instance, Mikhailov et al. 17 have used a co-recursion 
operator to generate an infinite hierarchy of cosymmetries for the Viallet equation, so a simple test to 
show that these produce an infinite hierarchy of CLaws would be useful. 
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